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Abstract 



We analyse dependence of the partition function on the boundary condition 
for the longitudinal component of the electric field strength in gauge field 
theories. In a physical gauge the Gauss law constraint may be resolved ex- 
plicitly expressing this component via an integral of the physical transversal 
variables. In particular, we study quantum electrodynamics with an exter- 
nal charge and SU(2) gluodynamics. We find that only a charge distribution 
slowly decreasing at spatial infinity can produce a nontrivial dependence in 
the Abelian theory. However, in gluodynamics for temperatures below some 
critical value the partition function acquires a delta-function like dependence 
on the boundary condition, which leads to colour confinement. 
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* * * 



It is with a feeling of deep sadness and irrecoverable loss that I had to com- 
plete writing this paper on my own. In the name of N.A. Sveshnikov Theoretical 
Physics has lost an extremely talented scientist, a great educator and a very nice 
person. His ideas and elegant scientific style had a big influence on his students, 
one of which I was privileged to be, and on many colleagues who knew him well. 
Perhaps, this paper may be viewed as a logical conclusion to years of our joint 
research on the problem of confinement. It was one of N.A. Sveshnikov's gifts 
to be able to see farther ahead in search for the hidden mathematical beauty of 
the physical world, the beauty that expresses the essence of all things. 



There is a considerable tradition in the field theory (see e.g. Ref. 0) to neglect any sur- 
face terms inevitably appearing in derivations. This is usually motivated by a fast decrease 
of all fields at spatial infinity. Such behaviour is only natural for theories with short-range 
interactions, but it is by no means obvious if long-range interactions come into play. 

The same problem acquires a somewhat different form at finite temperature. Obviously, 
the partition function of a translationally invariant system is ill-defined in infinite volume. 
Therefore, initially one has to enclose the system into a finite domain V and to assume 
some boundary conditions at the boundary dV. Then, to analyse dependence on boundary 
conditions we have to consider the functional, 



where Hy is the Hamiltonian of the system in volume V, (3 is the inverse temperature, is 
some subset of the canonical variables and function x defined on the boundary dV specifies 
the Dirichlet boundary conditions for the latter. In this framework we shall call Z[x] the 
effective partition function, its introduction being analogous in spirit to that of the effective 
action. It is reasonable to expect that Z[x], and any thermodynamic function, become 
independent of the particular choice of boundary conditions in the thermodynamic limit 
V — > oo. Once again, one may doubt whether that would really be the case for long-range 
interacting systems. 

It has been early recognised in the theory of gravity || that there are important physical 
situations in which boundary terms may have a physical meaning. We would like to mention 
also that many problems in hydrodynamics, such as e. g. description of surface waves 0, do 
require to consider variables at the boundary and nonvanishing surface terms. 

Recently the interest to boundary effects in various field theories has been rapidly in- 
creased 0. It has been found in the framework of the algebraic quantum field theory || 
that certain 2-dimensional models possess a nontrivial dynamics of the variables at infinity, 
and that such dynamics is responsible for the phenomenon of the dynamic mass generation. 

Our current purpose is to emphasise the role of boundary terms in the 4-dimensional 
gauge field theory at finite temperature and to study physical effects they can produce. In 
this case the appearance of a nontrivial x-dependence in Z[x] in the thermodynamic limit is 
almost obvious for the following reasons. The gauge theory is initially formulated in terms 
of an enlarged set of variables, the vector potentials A and the electric field strengths E in 
the Hamiltonian formulation, that make the gauge invariance explicitly manifest. Further, 
to obtain observable quantities one has to project the theory onto a subset of the physical 
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variables by resolving the Gauss law constraint, 'VE = p, and by adopting a gauge condition. 
Boundary conditions, of course, have to be compatible with these. For instance, in the 
Abelian case consider the boundary condition for the electric field strength component En 
normal to the boundary dV, which we take for simplicity as a sphere of radius R, 

R 2 E ll (Rx)\ dv = X (x). (2) 

This variable then has to obey the integrated form of the Gauss law, and therefore, 

Z[x}(x5( dxx- dxp). (3) 
JdV Jv 

We would like to emphasise that the dependence of the effective partition function on the 
boundary condition imposed on E\\ is of primary importance in the gauge theory because 
of the direct relation of this component to the colour charge flux due to the Gauss law. 
Analysis of namely this dependence will be the subject of the current paper. Based on the 
knowledge of this dependence alone we can suggest a simple confinement criterion: 

Z[ X ]^Y[S(x(x)). (4) 

X 

The latter condition simply means that the colour flux is strictly zero in every direction at 
spatial infinity for any state belonging to the Hilbert space of the system. 

The plan of the paper is as follows. In Sec. 2 we proceed with a careful calculation of 
of Z[x] in the simplest case of the Abelian theory with an external charge density. Then, in 
Sec. 3 we reproduce the answer obtained in the previous section using another technique, 
which is also applicable to the non- Abelian theory. Sec. 4 is devoted to calculation of Z[x] 
using the mean-field approximation and to consequent analysis of the confinement phase 
transition in gluodynamics. 



II. PARTITION FUNCTION OF QED WITH AN EXTERNAL CHARGE 

The partition function in a finite domain V may be represented by an Euclidean path 
integral over the fields periodic in time on the interval [0,/3], where (3 is the inverse temper- 
ature. The path integral is well-defined only if some boundary conditions are specified on 
all fields at the boundary dV. 

Let us denote the transversal and longitudinal components of vectors with respect to 
the gradient d as superscripts dA L = 0, and with respect to the vector x as subscripts 
x Aj_ = 0. For simplicity we choose the Coulomb gauge, A = A . We note also that the 
condition A = A± corresponds to the Fock-Schwinger gauge || (see Appendix A). We shall 
also assume that V = {x : \x\ = R} is a spherical domain with the radius R. The Gauss 
law constraint, 

8E = p(x) (5) 

allows to eliminate one space component of the electric field strength E. 

In a previous work |7j] on the basis of the general result || we have developed the Hamil- 
tonian formalism for the system in a finite spherical domain incorporating the boundary 
values as Hamiltonian variables. We have shown that the boundary conditions Ej_(Rx) = 
and XjFij(Rx) = are consistent with the localised time evolution in the Fock-Schwinger 
gauge. By transforming the theory to the Coulomb gauge one would arrive instead at the 
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boundary condition of the form E ± (Rx) = 0. Since these variables are independent, the 
dependence on a particular choice of the boundary conditions disappears for infinite system. 
The situation is quite different for the component E\\(Rx) = xE(Rx). Indeed, Eq. ([5]) may 
be easily solved, 

E ll (x) = \ [ X y 2 dy(p-dE ± )(yx). (6) 
x z Jo 

Both types of the transversal variables are connected by the relation, |J 

E ± (x) = E 1 - (x) -d r dy xE 1 - (yx) . (7) 

Jo 

Combination of Eqs. (H) and (^) now yields, 

R 2 E\\(Rx)- [ R y 2 dyp(yx) = A I* \r - y)dy xE L (yx) , (8) 
Jo Jo 

where we have used the spherical part of the Laplacian A = x~ l d 2 x + x~ 2 A. It is clear 
that the requirement E\\(Rx) = X (x), where xi x ) is arbitrary, is nothing but a constraint 
on the physical variables E x . As we have seen [7), it is this constraint that makes the finite 
volume Hamiltonian formalism closed. 

Now then, the partition function of QED with an external charge may be represented by 
the following path integral in the Coulomb gauge, 



Z = J VAVE 8(8 A) 5(R 2 E ll (Rx) - X (x)) 



exp J a d 4 x (i EA - l -E 2 + i AAA - ^(dE - p) 2 ) , (9) 

where we have used the notation for the domain A = [0, /3] x V (we shall also use the 
notation dA = [0,(3] x dV). We have also introduced a regularisation of the Gauss law by 
the parameter e, which should be set equal to zero at the end of calculations. This Gaussian 
integral is evaluated in a standard manner by a shift of the integration variables. To find 
the integral over E we introduce a new integration variable, Ei, 

E = E l + S, E 1 (Rx) = 0. (10) 

Here the new variable E\ satisfies the zero boundary condition and S is chosen so that there 
is no linear term in E\. This gives the equation on £, 

iA-£ + \d(dS)- \dp = 0, (11) 
R 2 xS(Rx) = X (x) . (12) 

The latter boundary condition follows from the second delta-function in Eq. (Q). We may 
decompose this vector onto the transversal and longitudinal parts in the momentum space, 
S = £ L — dip. Then, the transversal part is simply £ ± = iA, and the equation for <p 
becomes, 

(A-e 2 )^ = -p, (13) 

r2 m = ' X ^ ] ■ (14) 
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The partition function (^) is further decomposed as the product, 
Z = Z X Z , 

Z X = J VA^VE exp jf d'x ^-A 2 + ±AAA - l -E\ - ^(^i) 2 ) , (15) 
Z = expp(^-f dxx{x)<p(Rx) - - f dx p(x)^(x)] , (16) 

V £ J 8V Z JdV / 



where ip is the solution of Eqs. (0,0)- 

The solution of Eqs. fll3l , |14l) is, obviously, the sum of the homogeneous part, 0, satis- 
fying the nontrivial boundary condition, and of the inhomogeneous part satisfying the zero 
boundary condition, 

v = <P-G»p G=(A-e 2 )" 1 , (17) 

where G is the Green function corresponding to the zero Neumann boundary condition at 
R. The effective partition function can be presented as, 

Z = Z x Zp P Z px , 

Z x = exp(^ f dxxixMx)) , (18) 



2 JdV 

Z pp = exp J^dxdyp(x)G(x,y)p(y)j , (19) 
Z PX = exp (-^ J^dx<f>(x)p(x)j exp (~J^dxx(x) (G«p)(a;)^ . (20) 

It is natural to consider the problem further in terms of the spherical coordinates. Solution 
of Eq. flT3p regular inside the sphere is given by, 



Cu 



V^w«o> ( 21 ) 



v 'ler ■ - ■ 

n 




cosh^. (22) 



The constant C/ m is determined from Eq. (fl^). The role of the regulator e now becomes 
clear. The zero mode solution is simply, 



sinher _ xoo (c? ^ 



and it is 1/e 2 singular as e tends to zero. At the same time, solutions for other modes are 
perfectly regular in this limit and tend to, 

<f>im = C lm r l , Clm = -J R lk- ( 24 ) 

In consideration of the zero mode one must therefore be more careful and keep e nonvanish- 
ing. The zero mode Green function, that is defined by, 
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-Tr^r-e G 00 r,r = — , 25 



dG m (r,r') 



dr 

is easily calculated, 

I /I 1 



= 0, (26) 

r=R 



sinhei? — coshe-R/ei? . A 

H : — r-: — - , - smh er smh er . (27) 

cosh eit — smh eR/eR ) 



The leading terms at small e are, 



„ t a 3 9 lr 2 + r /2 1 

GooCr,/)^-— 3 + — --— -3 ? -, 28 

e^it' 3 5it 2 R 6 max(r,r) 

^3 3 r 2 \ 

+ 7^ Xoo • (29) 



e 2 i? 3 lOi? 2i? 3 



Substitution of these results into formulae fll8| - P0|) gives for the nonzero modes in the limit 

£ = 0, 

^>o = exp(-A ( £^!). (30) 

For simplicity, we assume that the distribution of the charge p is spherically symmetrical. 
In this case only the zero mode term survives. If we introduce the charge density momenta, 

rR rR 

Qoo = r 2 drp 00 (r), Goo = r A dr p 00 (r) , (31) 
Jo Jo 

our results may be summarised as follows, 

v a( 1 f R 2 j 2> Poo(x)poo(y) 3 2 

Zoo=eXP/3 l-2io XdXydy max(x,y) ~ 2^ (g °° ~ *° o) 
19 2 1 2 3 1 \ 

+ ^10^ 00 ~~ 10 ^ 00 _ Yq^ 00 ^ 00 ) ~ 2~R 3 "^ 00 ~ ^ 00 ^ 00 ) • ( 32 ) 

In the limit e — > this functional contains the delta-function of the condition Qoo — Xoo, 
and in addition we find the following correction to the standard answer due to the surface 
terms, 

Z 00 = exp (^jftQooj S (Qoo ~ Xoo) • (33) 

As a simple illustration of the above result it is instructive to consider the charge density 
poo = ~ K /{V^ r ) corresponding to a linear rising electric potential (p — nr. For such an 
exotic charge distribution, modeling the confinement alike potential, we find the additional 
constant contribution to the free energy density, 

It is interesting to note that this correction makes the free energy density smaller, and in 
this sense the boundary effects are thermodynamically significant. This example exhibits a 
promising connection between the boundary effects and the confinement phenomenon. 
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III. COLLECTIVE VARIABLE FORMULATION 



In the present section we give a different formulation of the same problem by introducing 
the collective variable, a, conjugate to the Gauss law constraint. Both formulations are 
completely equivalent in the Abelian theory, the transformation between them being just a 
trivial change of variables. However, the new formulation appears to be more fruitful in the 
non-Abelian gauge theory. Let us rewrite formula (|9|) in terms of the collective variable a, 
introduced by the definition, 



exp (-^ J^x(dE - p) 2 ^j = Jva exp J^d 4 x (^~a 2 + ta(dE - p)j. 



(35) 



One starts by taking the integral over E, 

J = J VE exp J^x (~^E 2 + iEA + iadE^j 5(R 2 E\\(Rx) -*(£)). (36) 

This can be done in analogy with the previous section by applying a shift E = E\ + £ and 
using the decomposition £ = S L — dip. Then, we shall get: £ ± = iA and ip = ia. Using the 
delta-function in the above formula one may rewrite the boundary term as i J dA dt dx x&. In 
the term / dx Ada the appropriate boundary term vanishes by the gauge condition. Thus, 
we obtain, 

J = exp jT d 4 x(A 2 + (da) 2 ) + i J^dt dx xcr^j 6 ( tR2 ^ + x)- (37) 
At this stage the \ dependence is contained only in the integral, 

Z = [Va exp f - f d A x {\{daf + %-a 2 + tap) 
J \ J a 2 2 

+ *J dA dt dx X crj ^R 2 -^ + \) ■ (38) 

Such integral is evaluated straightforwardly by a shift of the integration variable, a = a\+q. 
The new variable satisfies the trivial boundary condition, R 2 dai(Rx)/dR = 0, and q is 
found by the requirement that there is no term linear in o\ after the shift. This gives the 
equations on q, 

(A-e 2 )q = l p, R 2 ^ = *X- (39) 
Using the boundary condition on q we finally derive, 

Z = e W (- % -(5 J v dxqp+ % -(5 J dxqx) ■ (40) 

Note that with the identification q = —iip we have precisely reproduced the answer of Eqs. 
(0,0,0). 
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IV. NON ABELIAN THEORY 

In a previous paper || we have derived a representation for the partition function Z 
of SU(N) gluodynamics in a finite volume as the path integral over the collective variables 
analogous to that of the previous section. In order to find the effective partition function 
Z[x] dependence we have to return to the beginning of that derivation in the Fock-Schwinger 
gauge, 

Z[ X ] = f VAVEVa expf f d^x{iEA - -E 2 - -B 2 
J \J a 2 2 

+ iaVE)^ 5(A\\)5(R 2 E\\{Rx) - x ), (41) 

where a, clearly, is just a different notation for the temporal component of the gauge field, 
Aq. Obviously, any dependence on x is concentrated in the path integral over E\\, 

1 = JvEwexp^J^xi-^E 2 +ia(dx)E ll ))5(R 2 E ll - x ). (42) 

This is calculated by a shift E\\ = Ej + S, where S = —ida/dx and El satisfies the zero 
Dirichlet boundary condition. Such a derivation gives, 



(fP T 1 f f U(j\ f ~\\ OC 

L dt [-2Jv dx {T x ) + i J^x\)S(R 2 M -ix). (43) 

Next, by introducing the path integral representation, 

6XP {~lL dixB ^) = J Vu ex v(J A d4x (~]i u2 + iuB \\)) > ( 44 ) 

and after performing the integrations over A± and E± (see Ref. for more details) we 
finally obtain, 

Z[ x ] = [voVv exp(-W[<7, v\ + i [ dt dx ax) 5(R 2 o' - ix) , 

J J OA 

2W[a, v] = v • v + da • do + K_ • C+ 1 • K + 

+ K + • CZ 1 •K_ + ti log , (45) 

C± = —A x -V 2 ±D, K± = 8±v ± V t d±o , 

Vf = 5 ab d t - gt abc a c , D ab = gt abc v c , 

where the projected derivatives are defined by, 

di = rrj? dj , itl = hs ij - x { x j ± ie l]k x k ) (46) 

and the bullet denotes the 4-d integration over the domain A. 

In the saddle point approximation one expands the action near the saddle point, 

r a 5W 
W[q + ax] = W[q\ + / d 4 x —— ai (x) 
Ja oq{x) 

+ / dtR 2 dxS il) k]a 1 (Rx) + .... (47) 

JdA 



In the zero (mean-field) approximation we may write, 



Z\x] = exp (-Wfc] +ij d dtdx q X ) , (48) 

where the contribution from the first Euler derivative of the action £W [q] is precisely canceled 
with that from the surface term in Eq. (45). 

Thus, in the mean-field approximation the dependence Z[x] is controlled by the saddle 
point solution q. As we have seen in the previous section, the Abelian gauge theory possesses 
only the trivial solution q — 0. 

In Ref. H we have studied constant solutions of the mean-field equations. Let us re- 
produce those results briefly here, but in addition carefully keeping a finite volume. For 
simplicity we also restrict ourselves to the gauge group SU(2). First of all, we can introduce 
the notations for the free energy density, T R , 

W r = (3V r Fr, Fr = 1rF r , 1r = 8 ^ (0) , (50) 

P VR 

where V R = 47nR 3 /3 is the domain volume and 6(0) = 4- Z)j(2/ + 1) is the ultravioletly diver- 
gent angular delta-function with coinciding arguments. The function F R is now expressed 
via the dimensionless variables, 

2-k s . (2ixu\ 2 

pg \pg 



where to produce a real mean magnetic field v has to be purely imaginary (see Eq. 
After introducing the control parameter a = (27r) 4 / ' (2g 2 f3 Ar ) R ) and carrying out some deriva- 
tions we obtain, 

F R [u, s] = -au A +U R [u, s) , (52) 

U R [u,s]=U R [s}+V R [u,s}, (53) 

v r , , L R ((n + s) 2 + u 2 ) L R ((n + s) 2 — u 2 ) 

MuA = _ Ejog j^-^ , (54) 

,, , , P 4^ , {, COS27TS \ 

= ^S'° g i 1 " cosh( T (m + l/2), 8 /fi) J ' (55) 

where L R (x) = cosh(27iRy/x/ (3). 

It can be seen that at finite R this free energy possesses only a trivial minimum at s = 
u = 0. The situation changes after taking the thermodynamic limit, R — > oo. The resulting 
expression for the free energy density (see Ref. 0) possesses only the trivial stable solution 
u = s = at high temperatures. However, at some critical temperature T c the system 
undergoes a first order phase transition, below which there appears a deeper nontrivial 
minimum at u = s = 1/2 (see Ref. |§). 

Therefore, at high temperatures T > T c the dependence Z[x] is determined by the 
solution of the linearised equation 5 2 W/5q 2 • 0\ around q = 0. This can only produce the 
dependence akin to the Abelian theory. Namely, it contains the delta-function expressing 
the conservation of the global charge (Eq. @ with p = 0), and apart from that it is trivial in 
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the sense that Z[xi m ] — > 1 in the thermodynamic limit R — > oo. This situation, obviously, 
corresponds to the deconfinement phase as there is no restriction on the colour fluxes at 
infinity. 

However, below the critical temperature T < T c there is a nonzero constant solution 
kl = Tf/ofl- Since the system in invariant under the group of the big gauge transformations 
Goq parametrised by matrices U(x), the corresponding unit colour vector s(x) is arbitrary 
in every direction x. After integration over the orbits of the group SU(2) at each cone x 
the dependence becomes, 

z [x\ = [[— — 7^ ex P r dxx(x)\ 56 

where we have introduced a discretisation of the unit sphere with A being the infinitesimal 
cone area. It is well known |l[ that in the continuous limit the bare coupling constant 
vanishes g — > thereby making the effective partition function Z[x] a very sharply peaked 
function around the zero argument due to its essentially non-perturbative dependence on 
g. This property corresponds to the confinement phase, in which colour fluxes are equal to 
zero in every spatial direction at infinity. 

So, we can conclude that the dependence of the effective free energy on x is the following, 

z[x] = ^U & S(x(x)),T<%, (57) 

Therefore, the Gibbs average of an observable A is given by, 

<A> -\fd X (x)(A) x , T>T C , ^> 
where we have introduced the averages over inequivalent representations, 

(A) x = J T>aVpe- w[(7 ' v]+i IaA dtd& ^A[a,p] . (59) 

It is straightforward to see that any Gibbs average at low temperatures contains the singlet 
projector of the group G^, 

(A) = lim -LrrTr (e-? H » 8(Q R ) A) = lim -i-Tr (e~^ P s A) , (60) 

R^oo Zr[Ij\ R^oo Z/j[0J 

where Qr = J v dx V-E is the operator of the colour charge in volume V and P s is the 
singlet projector of the big gauge transformations. Presence of this projector in the Gibbs 
averages has been demonstrated to lead to the area law for the Wilson loop [H]|J, what is 
considered to be a standard confinement criterion. 



V. CONCLUSION 

The Gauss law in the gauge field theory may be resolved explicitly in a physical gauge. 
This produces effectively non-local interactions generating a boundary nontriviality of the 
theory. 

We have studied the dependence of the effective partition function on the Dirichlet 
boundary condition R 2 E\\ (Rx) = x( x ) imposed on the residual component of the electric 
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field for electrodynamics with an external charge and SU(2) gluodynamics. In the Abelian 
case this dependence always contains the delta function expressing the conservation of the 
total charge, and it is nontrivial only for charge distributions slowly decreasing at spatial 
infinity 

The non-Abelian self-interactions lead to a more unusual effect. Here the restriction of 
possible boundary values of the longitudinal electric field at low temperatures provides the 
confinement mechanism proposed by us in Ref. ||. Indeed, this quantity is proportional 
to the flux of the electric field through an infinitesimal cone in the direction x at infinity. 
Therefore, since the colour flux vanishes for any direction, no colour could escape to infinity 
and be experimentally observed. 
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APPENDIX A: GENERALISED FOCK SCHWINGER GAUGE 

Let Vr be a regular domain in R 3 topologically equivalent to a ball with a smooth 
boundary OVr. One can choose the curvilinear coordinate system X in Vr such that on the 
boundary OVr the first coordinate is constant and equal to some parameter R, which would 
play the role of an infrared regulator, i.e. OVr = {X : Xi(x) = R = const}. The field of 
vectors normal to the boundary for all possible values of R forms a differentiate vector field 
in R 3 . The ordinary Cartesian coordinates we shall denote as x to distinguish them from 
X . The local orthonormal frame then can be written as, 

where hj, are called the Lame coefficients. Also, to distinguish the components of vectors in 
the curvilinear frame we shall use the notations with indices in parenthesis, 

4» = 4»*. s «"ia^- (A2) 

Then em defines the field of normal vectors we have just introduced. It is natural to 
introduce the (2+1) decomposition onto the longitudinal and transversal (denoted by the 
Greek characters) components: i —>■ (l,a), a = 2,3, where the corresponding 2-d radius- 
vector will be denoted as X = (X 2 ,X 3 ). 

The gauge theory in a finite domain Vr, acquires especially elegant formulation in a gauge 
which is consistent with the shape of the boundary. Namely, we shall require that the normal 
component of the gauge field vanishes in every point, 

e (1) (x)A(t,x) = 0, A = A ± , A± = PA, P = 1 - e (1) (g)e (1) . (A3) 

This gauge condition is natural to name the generalised Fock-Schwinger gauge. Its most 
significant property is that the Gauss law, Vj-Ej = 0, can be resolved explicitly, 
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E(i) = ~\ /Jo) dX 'i ( h ®±)( X 'i, *)> ( A4 ) 

expressing the longitudinal component of the electric field through the transversal compo- 
nents of the gauge and strength fields, 

$± = ViE ±i = ^V (a) (J-E {a ^J . (A5) 

Above the lower integration limit xf^ is equal to some constant, which should be chosen in 
the reference point of our coordinate system. Analogously, one can solve the identity, 

XJ. Gi = , d = VjFij, (A6) 

expressing G(i) as, 

G(i) = ~\ Q } dX[ (h V t G ±t )(X[, X). (A7) 
Further, one of the components of the Bianchi identity is, 

Vj-Bj = 0, B k — -ei jk Fij, (A8) 
what allows us to express the longitudinal magnetic field as, 

= ~\ Q dX[ (h V t B ±l )(X[, X). (A9) 

Other components of the Bianchi identity e^fc e l / a \ V jE^ = give, 

1 d 



h a E {a) ) = V (a) (hiE (1) ) , (A10) 



A nice property of the Fock-Schwinger gauge is that the gauge strength and potentials 
are related to each other in a simple way. Indeed, from the definition of the gauge strength 
tensor applying the gauge condition we have, 

1 d 

Y a dxS KA(a) ^ = hlF(am ' (AU) 

However, the integral form of this equation, 

A (a) = ±- C dX[ (hihcFwXXl, X), (A12) 

flex J X-^ 

breaks the residual gauge transformations symmetry subgroup, whereas choosing the lower 
integration limits in all previous integral relations did not violate such a symmetry. Fixing 
some kind of the boundary condition above, e.g. choosing x[ ^ = x\ in the reference point, 
so that 
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\im (h a A (a) )(X) = 0, (A13) 

can be shown to be sufficient for determining a unique gauge field satisfying the gauge 
condition. Indeed, let A be any gauge field. The transformation to the Fock-Schwinger 
gauge is accomplished by a gauge transformation U(X), 

A(X) = U- 1 (A(X)-g- 1 d)U{X), E(X) = U~ 1 E{X) U(X), (A14) 
which can be found from the equation, 

1 d 



h x dX : 



V(X)=g(e {1) A)(X)U(X). (A15) 



For the moment we can choose the initial condition simply as U(X\ = 0,X) = 1. The 
solution of (| A 1 5|) is given by the Dyson P-exponent, 

U(X) = Pexp f 1 daR(a,X), R(a, X) = £/X 1 (/i 1 e ( i ) A)(aX 1 , X). (A16) 

J 

This can be explicitly worked out for the gauge fields, 

A\X) = A a (X) P exp f 1 da (-gt abc R c (a, X) ) - 

Jo 

-g~ l f 1 dpdR a (P,X)Pexp [ d-f {-gt ahc i? c ( 7 , X)). (A17) 
Jo Jo 

By applying an additional residual gauge transformation we can always satisfy the boundary 
condition (|A13| ), and uniquely. Really, suppose there exist two distinct gauge fields A' and 
A" satisfying Eqs. (TA"3| , |A"T3D and such that A" ^ U~ l A'U, V U = const. Then there should 
exist a gauge transformation between the two, 

A"(X) = U-\X) (A'(X) -g- l d) U(X) . (A18) 

Multiplication of this by e (1) gives (l/h 1 )(dU(X)/dX 1 ) = 0, i.e. U = U(X). Such 
transformations form the subgroup G res of the residual gauge symmetries in the Fock- 
Schwinger gauge. In fact, the boundary condition ( |A13|) does not permit such transformation 
since in the limit X — > x° in Eq. ( |A18| ) the gauge fields would have the singularity which 
is not compatible with such a boundary condition. This proves the uniqueness of the gauge 
field. 

As a simple example let us consider the elliptic coordinates: l<A A 1 <oo,— 1<X 2 <1, 
< X 3 < 2n, where X 3 = <fi is the polar angle and X\ = (ri + r 2 )/2a, X 2 = {r\ — r 2 )/2a. 
Such a coordinate system is defined by two points located at the distances ±a from the 
reference point along the z axis with r*i, r 2 being the radius-vectors from these points. The 
reference point is here x\ = 1, x\ — and the Lame coefficients are given by, 



hi = a 2 ^^, hi = a 2 ^—^, hi = a\Xl - 1)(1 - X|). (A19) 



In the limit a = the ellipsoid becomes a sphere and the gauge turns into the standard 
Fock-Schwinger gauge. For the spherical coordinates X 1 = r, X 2 = <fi, X 3 = 9 the Lame 
coefficients become very simple, 

h x = 1, h 2 = X x sin X 3 , h 3 = X v (A20) 

A technically attractive property of this particular gauge is that em do not depend on X\ 
and the vector normal to the boundary is just equal to the unit radius-vector em = x. 
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